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Abstract 

The disk fraction, the percentage of stars with disks in a young ciuster, is wideiy used to 
investigate the iifetime of the protopianetary disk, which can impose an important constraint 
on the pianet formation mechanism. The reiationship between the decay timescaie of the disk 
fraction and the mass dissipation timescaie of an individuai disk, however, remains unciear. 
Here we investigate the effect of the disk mass function (DMF) on the evoiution of the disk 
fraction. We show that the time variation in the disk fraction depends on the spread of the DMF 
and the detection threshoid of the disk. In generai, the disk fraction decreases more siowiy 
than the disk mass if a typicai initiai DMF and a detection threshoid are assumed. We find that, 
if the disk mass decreases exponentiaiiy the mass dissipation timescaie of the disk can be 
as short as IMyr even when the disk fraction decreases with the time constant of ~2.5Myr. 
The decay timescaie of the disk fraction can be an usefui parameter to investigate the disk 
iifetime, but the difference between the mass dissipation of an individuai disk and the decrease 
in the disk fraction shouid be properiy appreciated to estimate the timescaie of the disk mass 
dissipation. 
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1 Introduction 

Protostars are surrounded by disks consisting of gas and dust 
(referred to as a protopianetary disk), which have been investi¬ 
gated as sites of planet formation. Young stellar objects (YSOs) 
of low- and intermediate-mass, systems consisting of a protostar 
and a protopianetary disk, are roughly categorized into three 
classes by their spectral energy distributions (SEDs). Class I 
objects are the youngest with infrared emission from dust dom¬ 
inating the SED; Class II objects have intermediate ages with 
the SEDs characterized by a combination of stellar and hot dust 
emission; Class III objects are the oldest, consequently show¬ 
ing the weakest dust emission (e.g., Williams & Cieza 2011; 
Bachiller 1996). The mass and surface density of protopiane¬ 
tary disks decrease with time (Wyatt et al. 2007; Wyatt 2008), 
whose mechanism is still under debate. Although the dissipa- 
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tion involves a number of processes, photoevaporation likely 
plays an important role, especially for low-mass stars (Armitage 
2011). Planets are thought to form in the protopianetary disk. 
The dust and gas surface densities of the protopianetary disk 
are important parameters that determine the accretion rates to 
the cores and the final mass of gas giant planets to form (e.g., 
Ida & Lin 2004). It is therefore important to understand the evo¬ 
lution of the protopianetary disk for the investigation of planet 
formation. 

Different regions within a protopianetary disk can be inves¬ 
tigated with different wavelengths. Continuum emission in the 
near-infrared comes mostly from hot dust in the innermost part 
of the disk. Emission in the mid-infrared comes either from the 
surface of the disk or the mid-plane at the disk at a distance of 
several AUs. The far-infrared emission is dominated either by 
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the mid-plane or the outermost part of the disk (Dullemond & 
Monnier 2010). Atomic and molecular emission lines trace gas 
components and accretion activities. The investigation of the 
planet-forming disk evolution requires observations in the near¬ 
er mid-infrared. 

There are, however, difficulties in the observational investi¬ 
gation of the evolution of protoplanetary disks. First, the age 
of individual YSOs is difficult to estimate. It can be estimated 
from the loci in the Hertzsprung-Russel diagram (e.g., Strom 
et al. 1989; Kenyon & Hartmann 1995) or by near-infrared 
and X-ray photometry (Getman et al. 2014). However, these 
have inherent non-negligible uncertainties. The estimate based 
on the equivalent width of absorption lines has a better accu¬ 
racy (Takagi et al. 2010, 2011), but it requires high-resolution 
spectroscopy, which is not always available. The mass of pro¬ 
toplanetary disks is also difficult to estimate from near- and 
mid-infrared observations since they are optically thick at these 
wavelengths. 

The disk fraction, defined as the fraction of stars with disks 
in a star cluster, is thus widely used to study the evolution of 
protoplanetary disks since it relies on the age of star clusters, 
which reduces the uncertainties in the ages of individual stars, 
and since the presence of a disk can be estimated directly from 
the loci in the color-color diagram of near- and mid-infrared 
photometry as excess emission, although it does not estimate 
the disk mass. Haisch et al. (2001) show that the disk fraction 
estimated in the near-infrared decreases gradually with the clus¬ 
ter age, while Mamajek (2009) reports that the disk fraction de¬ 
creases exponentially with the time constant of 2.5Myr (decay 
timescale of the disk fraction, hereafter DTDF). 

The DTDF may change with the stellar mass and metal- 
licity of the cluster (e.g., Yasui et al. 2009). Hernandez 
et al. (2005) derive the disk fraction of nearby OB associa¬ 
tions using the JITFf-bands and suggest that the disk frac¬ 
tion for intermediate-mass stars is lower than that for low-mass 
(<1M0) stars. The result indicates the fast dissipation of the 
inner-disk for intermediate-mass stars. A similar result is re¬ 
ported by Carpenter et al. (2006), based on the Spitzer observa¬ 
tions of U Sco OB association. Kennedy & Kenyon (2009) in¬ 
vestigate the disk fraction of nine clusters, using infrared excess 
estimated from Spitzer observations and Ha equivalent width, 
and show that the disk fraction decreases with increasing stel¬ 
lar mass. Hernandez et al. (2008), using 2MASS and Spitzer 
data, find that the disk fraction of 7 Vel cluster is smaller than 
other clusters at a similar age (~ 5 Myr). They suggest that 
the small disk fraction may be attributed to the strong radiation 
field from massive stars in the cluster. For low-mass (< 2 M©) 
stars, the disk fractions estimated at the near- and mid-infrared 
do not show any appreciable difference between each other. On 
the other hand, Yasui et al. (2014) suggest that the disk dis¬ 
appears by ~ 3Myr faster at the near-infrared than at the mid- 


infrared for intermediate-mass (~1.5-7M0) stars, being con¬ 
sistent with the lack of planets discovered in the vicinity of the 
intermediate-mass stars. Ribas et al. (2014) investigate the disk 
fraction of 22 young clusters using SED fitting from optical to 
mid-infrared. They show that the primordial disk detected at 
shorter (3.4-12/im) wavelengths disappears faster than that at 
longer (22-24 pm) wavelengths, indicating that the disk closer 
to the star is evolving more rapidly. 

Although the evolution of the disk can be studied efficiently 
by the disk fraction, the DTDF is not equal to the mass dissipa¬ 
tion timescale of an individual disk (MDTID). The DTDF is a 
statistical measure and does not directly indicate the variation in 
the disk mass of individual objects. It is the timescale of the disk 
dissipation that governs the planet formation process. However, 
the reliability of the DTDF as an estimator of the MDTID has 
not been studied systematically. It is thus of importance to un¬ 
derstand the relation between the DTDF and the MDTID in a 
semi-quantitative manner. 

How the disk fraction decreases with the cluster age depends 
on the statistical properties of the cluster. Owen et al. (2011) 
investigate the photo-evaporation mechanism of the protoplan¬ 
etary disk by X-rays with radiation-hydrodynamic simulations. 
They demonstrate that the observed disk fraction at the near- 
infrared can be reproduced by the X-ray luminosity function of 
the cluster. Therefore, the behavior of the disk fraction must 
be investigated taking account of the distribution of physical 
parameters of the cluster. Sub-millimeter observations indicate 
that the total dust mass in the disk is distributed over a wide 
range (e.g., Andrews & Williams 2005, 2007). Armitage et al. 
(2003) show that the observed spread of the disk lifetime in the 
Taurus cloud is consistent with theoretical models of disk evo¬ 
lution if a dispersion in the initial disk mass is assumed. The 
initial disk mass function (DMF) can thus affect the time vari¬ 
ation of the disk fraction in a statistical way, but this effect has 
not been investigated in detail yet. 

In this paper, we devise a simple model to estimate the ef¬ 
fect of the initial DMF on the time variation of the disk fraction. 
The difference between DTDF and MDTID is calculated, tak¬ 
ing account of the broadness of the DMF. Details of the model 
are described in Section 2. In Section 3, the variation in the 
calculated disk fraction is investigated and compared with the 
observed disk fraction. The results are summarized in Section 4. 

2 Model 

The evolution of the disk fraction is investigated by taking into 
account the disk mass distribution. In the present model, we 
consider the observation of a young star cluster in the near- or 
mid-infrared, where the protoplanetary disk is optically thick. 
In observations, the target stars are selected by their spectral 
types and thus the stellar mass is in a narrow range. In the 
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following, we do not consider the distribution of the mass of 
the central star. Infrared excess (/®’‘) is assigned to each star. 
The star is recognized as being associated with a protoplanetary 
disk when is larger than a critical value (/**'). All stars in 
the cluster are assumed to be bom as a single star at the same 
time with different disk masses. The effect of the age dispersion 
within the cluster is discussed in Section 3.3. The disk fraction 
is defined by the fraction of such stars with > /*'' in the 
cluster. 

The DMF, the fractional number of the disk with the mass 
between ma and ma+dma, is defined by cl>o{'md)dm and 
4>o{rad) is normalized as 

J 0o(ma)dma = 1. (1) 

Based on a survival analysis, Andrews & Williams (2005) show 
that the DMF of the Taurus-Auriga region is well approximated 
by a log-normal distribution. Similar results are also reported 
by Andrews & Williams (2007) and Mann & Williams (2010). 
According to their results, the initial DMF is approximated by 


(?l)o(ma)dma = 
1 


ma 




- exp 


1 / logma-logM 

2 V a 


2 ' 

dma, 


( 2 ) 


where p and a define the location and width of the distribution. 

The disk dissipates with time. Luhman et al. (2010) suggest 
that the inner optically thick disk evolves into an optically-thin 
phase rapidly because the number of transitional disks may be 
small. The evolution of the inner disk mass is, however, not 
well understood. If the disk is regarded as an isolated system, 
the mass loss rate of the disk should be determined by the pa¬ 
rameters of the disk itself, such as the disk mass, the luminosity 
of the central star, and the angular momentum. To simplify the 
case, we assume that the dissipation of the disk mass is given 
by an unary function of the disk mass; 


dmd 

df 


(3) 


C("id)’ 

where C{m) is an arbitrary positive definite function to define 
the disk mass loss rate. By integrating Equation (3), the rela¬ 
tionship between ma (t) and t is obtained. 


r™d(t) 


(■(m)dm = —t. 


/rridfO) 


(4) 


We define Z (m) as the integration of ( (m). 


2 (md(f)) — 2 (md(0)) = —t. 


(5) 


Since C{m) is positive definite, the integrated function Z{m) 
is a monotonically increasing function, for which the inverse 
function Z~^{m) is uniquely defined. Thus md{t) is solved as 

md{t) = M (t,md(0)) = Z~^ (-t + Z (md(0))). (6) 

The inverse function Z~^(m) is also a monotonically increas¬ 
ing function. A4 (t,md(0)) monotonically decreases with t: 


t> t' <;=^ M (t,md(0)) < M (f',md(0)) . (7) 

The MDTID is defined as a typical timescale of the de¬ 
crease in A4 (t,md(0)) as discussed in Section 3.1. Instead, 
A4 (t, md(0)) is an increasing function in terms of md(0): 

md(0) > md(0) 4=^ M (t,md(0)) > M (f,md(0)) . (8) 

Equation (8) assures that the mass dissipation process does not 
change the order of the disk mass. Equation (5) can be solved 
in terms of md(0): 

md(0) = Z~^ {t + Z (md(f))) = M {-t,md{t)). (9) 

The relationship between the initial DMF and the DMF at t is 
given by 

(fit (md)dmd = (j)o {M (-f,md))dAl {-t,md) ■ (10) 

The DMF changes in time depending on the time evolution of 
M{—t,md). Although we assume that the initial DMF is given 
by a log-normal distribution, the shape of the DMF at t can be 
different from a log-normal function. Here we discuss a general 
case of M{—t,md). 

The near-infrared excess is directly related to the inner disk 
mass, not the total disk mass. The mass function of the in¬ 
ner disk is, however, not observationally constrained. Here, 
we simply assume that the inner disk mass is proportional to 
the total disk mass. The amount of the excess should depend 
on the disk inclination and the shape of the inner rim of the 
disk (Dullemond & Monnier 2010). The connection between 
and the disk mass has not yet been understood well. For 
sake of simplicity, we assume that is given by an unary 
function of m-d- The infrared excess /®’‘(md) should be a non¬ 
decreasing function with nid- The threshold mass is de¬ 
fined by f^^imfd) — The disk is detected in the infrared 
when nid > rnf^, while infrared excess is not detectable when 
rrid < In this formalism, we do not specify the relation¬ 
ship between and rrid- The presence or absence of the disk 
is determined only by rrid, irrespective of the functional form of 
/™(md). The infrared excess may depend on several pa¬ 
rameters such as the disk inclination. Justification of this simple 
formulation is discussed in Section 3.3. 

The disk fraction at the age of t is calculated by 

poo 

iF{t)= / 0t(m)dm. (11) 

JrrJ^ 

Then, by substituting Equations (2), (3), and (10) in 
Equation (11), the disk fraction T{t) becomes 

logTW (-i,md'’) -log/r] 

71^ J ’ 

where erfc(a;) is the complimentary error function. Since 
erfc(x) is a decreasing function on x, the disk fraction, 
if (t, cr), decreases as t increases. The DTDF is defined 

as a typical timescale of the decrease in Equation (12). 


T {t-,rrid,p,a) = -erfc 
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11 = 1x10 ^Msun, CT = 1.31dex 



Fig.l. Time evolution of the disk fraction against Note that the horizontal 
axis is not the time t, because the variation of the disk fraction can be 
described in terms of J in a more general manner than of t (see text). The 
disk fractions normalized at age zero are plotted for 
mj** = 2.Ox 10“"', l.Ox 10“"', and 0.5x 10“"' Mq by the solid, dashed, 
and dot-dashed lines, respectively. The median and variance of the DMF 
are fixed at 10“^ Mq and l. 3 ldex, respectively. The dotted curve shows 
the reciprocal of indicating the decrease in the individual disk mass (see, 
text). 


3 Discussion 

3.1 Behavior of Disk Fraction 

Since the complimentary error function erfc(2;) is a monoton- 
ically decreasing function, the disk fraction will decrease with 
increasing A4 (—t,md(0)), which increases monotonically with 
t. Equation (12) shows that the central loci of the DMF p shifts 
the disk function along the axis of age. The width of the DMF 
cr normalizes the argument of the complimentary error function 
in Equation (12). As rr increases by a factor of x, the time vari¬ 
ation of the disk fraction is decelerated by a factor of x. This 
suggests that the broadness of the DMF has a large impact on 
the evolution of the disk fraction. 

To investigate a typical decreasing timescale of Equation (6), 
we define the time T(^) such that the disk with the initial mass 
of loses its mass down to mf' at f = T{^): 

(13) 

The duration means an observational lifetime of the disk with 
the initial mass of Equation (13) is equivalent to 

^ no = 2 (Cm*") - 2 (mf) . (14) 

T{0 increases with ^ since Z{m) is an increasing function. 
When ^ is unity, T{0 is equal to zero. The time T{0 is re¬ 
garded as a typical timescale for the disk mass to decrease by 
a factor of 0 By substituting Equation (14), the following rela¬ 
tionship is obtained: 

M{-T{0,mf)=Z-^{n0 + Z(.^f))=^rnf. (15) 

The disk fraction at T{0 is obtained as 


T(T{0-,mf,p,a) = ^erfc 

The disk fraction decreases with ^ as well as T{0- 
Equation (16) does not include T{0 explicitly, indicating that 
the relationship between the disk fraction and ^ does not de¬ 
pend on the functional form of T{0- Thus, the evolution of the 
disk fraction can be described in a general manner in terms of 
0 Figure 1 shows the evolution of the disk fractions normal¬ 
ized at age zero against The DMF of the Taurus cloud is 
assumed to be a representative one, so that p and a are fixed 
af l.OxlO“®M0 and 1.31 dex, respecfively (the “full” sample 
in Andrews & Williams 2005) ^ By comparing sub-mm and 
near-infrared observations, Andrews & Williams (2005) suggest 
that the near-infrared (Ks—L) color excess method can detect 
the protoplanetary disk with a mass down to ~ 10“"^ Mq. The 
disk fractions are plotted for = 0.5x10“'*, 1.0 x 10“'*, and 
2.0x 10“* Mq. The gray dotted curve shows the reciprocal of 
0 representing the dissipation of the disk mass. If the disk frac¬ 
tion decreases with the same timescale as the disk mass, the 
curve should follow the dotted curve. When the disk fraction 
decreases slower than the disk mass, the curve will deviate up¬ 
ward, and vice versa. All the disk fractions deviate upward in 
Figure 1. An e-folding time of the individual disk mass dissi¬ 
pation estimated by T(e), which corresponds to the MDTID. 
On the other hand, the disk fraction for = l.Ox 10“* does 
not decrease by a factor of e until ^ increases to about 40, sug¬ 
gesting that the DTDF is as long as T(40). We found that 
the disk fraction decreases as fast as the disk mass only when 
mn ~ 2x10“* Mq. Such an extreme case is excluded from 
observations. Figure 1 indicates that the disk fraction generally 
decreases slower than the disk mass when a typical initial DMF 
and a detection threshold are assumed. 


log (C»TdVp) 


(16) 


3.2 Case Study: Exponential Decay 

To quantitatively compare the MDTID and DTDF, we consider 
a special case of C(m) = t jm, which corresponds to an expo¬ 
nential decay. Equation (6) is reduced to 

md(f) = md(0)e“^. (17) 

Here r is equivalent to T(e). In this section, we define the 
MDTID by the e-folding time r. With the assumption of 
C(m) = r/m, the DMF at t is given by 


o (md)dmd = 


1 

nid y/na 


exp 


_ t 

logmd—logpe T 


2 ' 


dmd. 


(18) 


' Although Andrews & Williams (2005) have reported that the median and 
the variance are 5x 1 O“'*M 0 and O.SOdexforthe “detection" sub-sample, 
we adopt the values for the "full” sample because non-detected samples 
should be properly taken into account. 
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Cluster Age 



Fig. 2. a) Disk fractions for different o-’s. b) Disk fractions for different 
The horizontal axis is the cluster age t measured by r. The red 
dot-dashed line shows n exponential decay with a time constant of t. 


Here, the DMF remains a log-normal distribution. The median 

_ 

of the DMF at t is given by pe -r, indicating that the location 
of the DMF shifts with the time constant of r, while the scale 
parameter of the log-normal distribution, a, is constant in time. 
By applying Equation (17), the disk fraction is given by 

f/r-hlog(md*'//i)1 

y2a J ■ ^ ^ 

Figure 2 shows the time variation of the disk fractions for dif¬ 
ferent DMFs. The horizontal axis denotes the cluster age t mea¬ 
sured by r. The median mass of DMF is set to be 1.0 x10“^Mq 
(A ndrews & Williams 2005). Figure 2a shows the results for 
different DMFs, <j = 0.65, 1.31, and 2.62 dex, where mj)*' is 
fixed at 10“'* Mq. The figure indicates that the disk fraction 
decreases slower as the disk mass function becomes broader, 
because t is normalized by ar in Equation (19). The disk frac¬ 
tion is not unity even at age zero, simply because some disks do 
not have a large amount of initial mass enough to be detected. 
Figure 2b shows the disk fractions for different detection thresh¬ 
olds, mf = 0.5x10“*, 1.0x10“*, and 2.0x10“* M©. The 
disk fraction shifts to the left with increasing m^''. The param¬ 
eter affects the disk fraction at age zero, but does not 

change the slope of the disk fraction. 

The red dot-dashed line in Figure 2a shows an exponential 
decay curve with the time constant of r: exp(—f/r). All the 
disk fractions in Figure 2a decrease more slowly than the expo- 


T (t,mf,p,a) = -erfc 



Cluster Age (Myr) 


Fig. 3. Comparison with observed disk fractions. The gray filled circles with 
errors are observed disk fractions from Mamajek (2009). The blue solid line 
shows the disk fraction with the best-fit parameters. The red dot-dashed 
line shows an exponential decay with a time constant of 2.5 Myr. 


nential curve. The e-folding time of the disk fractions ranges 
between 2r and 4r, depending on This is consistent with 
the discussion in Section 3.1: the disk fraction decreases gener¬ 
ally more slowly than the disk dissipates. Figure 2 implies that 
the difference between the MDTID and DTDF is not negligible. 

As a demonstration, we calculate the difference between 
DTDF and MDTID using a real data set. Since the present 
model is highly simplified, this practice is to demonstrate the 
degree of effect of the initial DMF on the estimate of the 
MDTID. We do not intend to derive an accurate estimate of the 
MDTID. Figure 3 shows that the observed disk fractions of a 
number of star clusters collected in Mamajek (2009) and refer¬ 
ences therein. In the current model, the disk fraction is given 
by Equation (19). The location and width of the initial DMF, 
p and a, are fixed at 1.0x10“® M© and 1.31 dex, respectively 
(Andrews & Williams 2005). The other two parameters, r and 
m}^, are estimated by a maximum likelihood method. The like¬ 
lihood function is defined as 


£(t,7) oc n exp 

i 


1 / 

2 M 


2-1 


( 20 ) 


where fi and Afi are the observed disk fraction and its un¬ 
certainty, and ti is the cluster age. The parameters, r and 
are estimated by maximizing Equation (20). The er¬ 
rors in r and are estimated by parametric bootstrapping. 
The best fit parameters are r = 0.921 q' 28 Myr and mj)*' = 
0.4810 25 xl0“*M©. Even though the model is simple and 
ignores many effects, the detection threshold presented above 
is consistent with the value reported by Andrews & Williams 
(2005) (~10“* Mq). The blue line with the crosses in Figure 3 
shows the fitted curve. Mamajek (2009) suggests that the evo¬ 
lution of the observed disk fraction is well represented by an 
exponential decay with the time constant of 2.5Myr, which is 
shown by the red dot-dashed line. The decrease in the disk frac¬ 
tion is well approximated by exp(—f/T), where T — 2.5Myr, 
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suggesting that DTDF is about 2.5 Myr. The present result, how¬ 
ever, suggests that MDTID r is as short as 0.9Myr. 

3.3 Justification of Modei 

The discussion in Section 3.1 is applied to a general case, where 
various physical processes — such as binarity, variable X-ray 
flux, and photo-evaporation — can contribute to the dissipa¬ 
tion. It assumes that the initial DMF is given by a log-normal 
function, but does not have a constraint on the DMF at a given 
time. Although the present model is rather robust in this sense, 
it is simplified and based on several assumptions. Effects of the 
simplification on the results are discussed. 

In the present model, the stars in the cluster are assumed to 
be coeval. The star formation is, however, not instantaneous. 
Getman et al. (2014) investigate the age of individual YSOs in 
massive star forming regions and show that the spread of the 
stellar ages in a cluster is about 1-2 Myr, comparable with the 
DTDF (~ 2.5Myr). The age spread within a cluster should not 
be neglected in the discussion of the disk fraction. At a given 
epoch, YSOs formed earlier have lost more mass than those 
formed later, and vice versa. Thus, the apparent width of the 
DMF is broadened if the age distribution is taken into account 
in the present model. While the width of the DMFs has been 
observationally measured (e.g., Andrews & Williams 2005), the 
intrinsic width of the DMF and the broadening effect due to the 
age spread are not distinguished. The observationally estimated 
width of the DMF should include the effect of the age distri¬ 
bution. Since we employ the observationally estimated width 
of the DMF (Andrews & Williams 2005) in our discussion, the 
present results indirectly take account of the effect of the age 
distribution. 

The present model assumes that the mass of the inner disk 
is proportional to the total disk mass. This is not observation- 
ally confirmed, because the inner disk is optically thick in the 
infrared and the infrared excess does not reflect directly the 
amount of the inner disk mass. The present results, however, 
mainly depend on the width of the initial DMFs. The results are 
expected to be valid as long as the distribution of the inner disk 
mass is relatively as wide as the DMFs observed in the sub-mm 
wavelength. 

The present model does not include any inclination effect 
on the infrared excess, although the infrared excess should 
be changed with inclination. The relationship between 
and the inclination angle may be complicated. Since the near- 
infrared emission is mainly from the inner rim of the dusty 
disk (Dullemond et al. 2001; Natta et al. 2001; Muzerolle et al. 
2003), the effect of the inclination angle can heavily depend on 
the shape of the inner rim (Dullemond & Monnier 2010). If the 
inner rim has a round shape (Isella & Natta 2005), the effect of 
the inclination on can be small. Observationally the infrared 


excess does not show any strong correlation between the incli¬ 
nation (Dullemond & Monnier 2010), suggesting that the incli¬ 
nation effect on the present result is not significant. Although 
the inclination significantly changes the infrared excess /®’‘, it 
is hard to distinguish the contribution from the inclination and 
the disk mass unless the inclination angle is determined. The 
effect of the inclination can be practically taken in the present 
model by broadening the initial DMF. 

The relationship between the infrared excess / and the disk 
mass, which is important to connect the initial DMF with the 
disk fraction, remains to be understood. We assume that the in¬ 
frared excess / does not decrease as the disk mass increases. 
In the present model, the infrared excess is given by an 
unary function of the disk mass. The amount of the infrared ex¬ 
cess is closely related to the shape of the inner rim (Dullemond 
& Monnier 2010; Isella & Natta 2005). The shape of the rim 
depends on the pressure scale height at the inner rim 
which is proportional to 1/x/M*, where M* is the stellar mass 
(Isella & Natta 2005). Therefore, the neglect of the mass de¬ 
pendence can make a moderate effect on the result and should 
be taken into account in the next step. The effect on the present 
results due to this simplification is not expected to be severe as 
long as the mass range is sufficiently narrow. 

In Equation (3), we assume that the decrease in the disk 
mass is given by an unary function of md- This assumption can 
be rephrased as disks with the same disk mass have the same 
lifetime. The decreasing rate of the disk mass should depend 
on the radiation from the central star and neighborhood stars. 
The assumption can become invalid if the mass range is not 
sufficiently narrow or the radiation field signiflcanfly changes 
within the cluster. Consequently, the lifetime of the disk can be 
extended or shortened. The effects of the mass range and the 
radiation field of neighborhood stars can be qualitatively esti¬ 
mated by broadening the width of the initial DMF. As shown 
in Equation (12), the disk fraction decreases more slowly than 
the disk mass as the initial DMF becomes wider. Although the 
present model does not correctly take account of the effects of 
the mass range and the radiation field from neighborhood stars, 
those effects should not change the conclusion of the present 
result qualitatively. 

3.4 Implications of Results 

As shown in Figure 2, the apparent evolution of the disk frac¬ 
tion depends on the shape of the initial DMF. The disk frac¬ 
tions of clusters with different initial DMFs should follow dif¬ 
ferent curves. In most observational studies, the disk fractions 
are assumed to follow a single curve (e.g., Haisch et al. 2001; 
Yasui et al. 2009, 2014; Ribas et al. 2014). In other words, 
they assume that the initial DMFs are assumed to be the same. 
However, the DMF has so far been obtained for only a few clus- 
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ters, e.g., in the Taurus, Ophiuchus, and Orion regions (e.g., 
Williams & Cieza 2011). Those observations suggest that the 
DMF can be approximated fairly well by a single log-normal 
function, while Mann et al. (2015) suggest that the DMF of 
NGC2024 is possibly top-heavy in comparison with that of 
the Taurus cloud and Orion Nebula Clusters. If the DMF is 
top-heavy in the sample clusters, the DTDF will become even 
longer than the present estimate. Only by investigating the evo¬ 
lution of the disk fraction, it is not possible to disentangle the 
variation in the MDTID and the variation in the initial DMFs. 
Further investigation about the variation in the initial DMF is 
needed. 

Previous studies have been compared the disk fractions of 
different clusters and discussed the DTDF (e.g., Haisch et al. 
2001; Yasui et al. 2010). The present results suggest that the 
apparent evolution of the disk fraction depends on the thresh¬ 
old mass. The disk detected in the near-infrared is expected 
to be fully optically thick. The amount of excess in the 
near-infrared is insensitive to the disk mass in the optically thick 
phase (Wood et al. 2002). The choice of the detection thresh¬ 
old (/‘*') will not matter. In general, the transition from a thick 
to thin disk occurs in a small mass range at the near-infrared 
wavelength. The fraction of disks in the transition in a cluster 
is small (e.g., Skrutskie et al. 1990; Wolk & Walter 1996; Cieza 
et al. 2007). Thus, the dependence of the threshold mass mj)** 
is negligible. At wavelengths longer than 25 pm, the disk can 
be detected in the optically thin phase and the excess becomes 
sensitive to the disk mass (Wood et al. 2002). The difference in 
the detection threshold may have to be taken into account. 

4 Conclusion 

The disk fraction, which is the fraction of stars with disks in a 
young cluster, is widely used to observationally investigate the 
lifetime of the protoplanetary disk. The time evolution of the 
disk fraction should depend on the disk mass function (DMF) at 
age zero. We discuss a simple model to analytically investigate 
the relationship between the mass dissipation timescale of an 
individual disk (MDTID) and the decay timescale of the disk 
fraction (DTDF). 

In the present model, the evolution of the disk fraction can 
be described by the detection threshold of the disk mass (mjf), 
and the loci (p) and the dispersion (cr) of the initial DMF. The 
DMF of the Taurus cluster has per 10“® Mq and a ~ 1.31 dex 
(Andrews & Williams 2005). Using the near-infrared color ex¬ 
cess method, the disk with the mass of ~10“"' Mq can be de¬ 
tected (Andrews & Williams 2005). Assuming these parame¬ 
ters, the present model indicates that the disk fraction generally 
decreases more slowly than the disk mass, suggesting that the 
DTDF is longer than the MDTID. Given that the disk mass dis¬ 
sipates exponentially, the difference between the DTDF and the 


MDTID is suggested not to be negligible. Comparison with the 
observational data (Mamajek 2009) suggests that the MDTID 
is as small as 1 Myr even if the disk fraction exponentially de¬ 
creases with the time constant of 2.5Myr, corresponding to the 
DTDF. Although the present model is simple and primitive, we 
are confident that the present results are qualitatively valid. 

The present results suggest that the evolution of the disk 
fraction depends on the shape of the initial DMF. The variation 
in the DMF remains to be understood. Further observational 
investigation is needed. Although the apparent evolution of the 
disk fraction may depend on the detection threshold mass, the 
matching of the detection threshold does not matter unless the 
disk fraction is measured at a longer (>25 pm) wavelength. 

The study of the protoplanetary disk is being accelerated by 
the advent of ALMA, which can resolve the disk structure and 
precisely measure the mass of the disk. However, until the disk 
mass is accurately measured for a sufficiently large sample of 
disks in various environments, statistical approach is still im¬ 
portant for the study of the disk evolution. The disk fraction 
remains to be efficient and useful if the relationship between 
the DTDF and the MDTID is properly appreciated. 
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